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ABSTRACT 

Objectives: To identify the differential effects of 
patient, health service, temporal and geographic factors 
on length of stay (LOS) for chronic obstructive 
pulmonary disease (COPD)-related admissions. 
Design: We used stratified Cox proportional hazard 
model to evaluate the association between LOS and 
patient, health service, temporal and geographical factors. 
Setting: Patients resident in Blackpool, North West 
England, admitted to the local hospital with COPD. 
Participants: We used the Admitted Patient Care General 
Episode Commissioning Dataset for the period 1 April 
2005-31 March 2010. We analysed records of admission 
spells among patients resident in Blackpool aged 40 years 
or older admitted with a primary diagnosis of COPD. 
Results: There were 2410 admissions meeting the 
inclusion criteria over the period. These admissions were 
attributed to 1172 COPD patients, an average of 2.06 
admissions per patient. The median LOS was 6 days 
(95% CI 6 to 6) while the mean was 9.8 days (95% CI 9.1 
to 10.5). Patients were 22% more likely to be discharged 
earlier in 2009/201 0 compared with 2005/2006 (adjusted 
HR 1.22; p=0.01 00). LOS was associated with 
socioeconomic deprivation with those in the most 
deprived areas being 35% less likely to be discharged 
earlier compared with those from the least deprived areas 
(adjusted HR0.65; p=0.0010). 
Conclusions: LOS among COPD patients have reduced 
over the period of the study. Age, deprivation, Charlson 
index, specialty of admission and cause of exacerbations 
were independently associated with LOS. Though there 
were no significant associations between LOS and season 
of admission and distance from hospital, there were 
significant variations in LOS associated with these 
variables based on selected patient characteristics. 
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INTRODUCTION 

Chronic obstructive pulmonary disease 
(COPD) is projected to become the third 
most common cause of morbidity worldwide 
by 2020. 1 In England COPD prevalence is esti- 
mated at between 2% and 4%. 2 During the 
2009/2010 financial year COPD admissions 



ARTICLE SUMMARY 



Article focus 

The focus of this study is on: 

■ The length of hospital stay (LOS) among patients 
admitted with chronic obstructive pulmonary 
disease (COPD) in Blackpool, a spearhead 
unitary local authority in the North West of 
England admitted to the local hospital. 

■ The associations between LOS and patient, health 
service, temporal and geographical factors. 

Key messages 

■ LOS among COPD patients have reduced over 
the study period. 

■ COPD patient from the most deprived areas 
stayed longer in hospital compared with their 
counterpart from the most affluent areas. 

■ Effect of season on LOS varied by the level of 
deprivation. 

Strengths and limitations of this study 

■ Our study provided a unique insight into the 
experiences of COPD patients from a small geo- 
graphical area using a local hospital. 

■ Our study takes into account case-mix and inter- 
action effects among predictor variables. 

■ Incompleteness and inaccuracy are the main pro- 
blems associated with routine data we used. 

■ We could not assess the effects of other potential 
confounders on LOS because the dataset we used 
did not contain information on these variables. 

■ We used population level measure of deprivation 
(Index of Multiple Deprivation) that could pre- 
dispose to ecological fallacy. 



accounted for about 1.1% of in-patient fin- 
ished consultant episodes and 1.5% of all fin- 
ished consultant episodes bed days. 3 

A high proportion of the costs of man- 
aging COPD is attributed to exacerbations. 
This varies from 40% to 57% of the total 
direct costs and can be as high as 63% in 
severely affected patients. 4 Inpatient admis- 
sion is a major cost driver in COPD manage- 
ment'' h accounting for about 54% of the 
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direct cost associated with its management in the UK 
The cost of hospitalisation is also significantly influenced 
by the length of hospital stay. 8 

Though many studies on improving admission outcomes 
for COPD patients and improving efficiency tended 
to focus on reducing duration of in-patient care 9 10 no 
optimum length of stay (LOS) has been consistently 
described. International trends in COPD-related hospital- 
isation have shown that although the average LOS has 
decreased since 1972, admissions rates have increased in 
patients aged 45 years or older. 8 This has led some to 
suggest that efforts to cut down the duration of inpatient 
episodes may result in a 'revolving door' phenomenon. 11 

According to the English Department of Health 2010 
COPD strategy consultation document 12 the LOS for 
COPD-related admissions has reduced to a median of 
5 days in 2008 from a median of 6 days in 2003 but 
re-admission rates have been increasing. This is despite 
the fact that relatively few patients with exacerbations 
are admitted. 13 Taking cognizant of these patterns, the 
consultation document identified a key objective of 
reducing 'length of stay for people with COPD, where 
appropriate'. 12 

Our objective for this study is to identify the differen- 
tial effects of patient, health service, temporal and geo- 
graphical factors on LOS for COPD-related admissions. 
Understanding of these influences will help optimise 
care for COPD patients needing hospital admissions. 

METHODS 

Study setting and predictor variables 

Blackpool is a spearhead unitary local authority in the 
North West of England with an estimated population 
of 139 974. 14 The majority of population live in areas 
classified within the fifth-most deprived national depriv- 
ation quintile. The Blackpool Teaching Hospitals NHS 
Foundation Trust (BTH) is the main provider of 
inpatient hospital care for Blackpool residents. 

In tandem with national and international trends, COPD 
prevalence in Blackpool continues to increase. Estimates 
based on general practice Quality Outcome Framework 
(QOF) data showed that the prevalence in Blackpool was 
2.7% in 2009/2010 up from 2.3% in 2004/2005. The 
2009/2010 figure was significantly higher than the North 
West average of 2.1% and England average of 1.6%. 

The predictor variables considered in the study are 
listed in box 1. These were selected because of their 
well-established association with health outcomes such as 
life expectancy, quality of life, uptake of preventive ser- 
vices and treatment outcomes, 16 health service use 17 
and seasonal variations in mortality. 18 

Data 

We used the Admitted Patient Care General Episode 
Commissioning Dataset for the period 1 April 2005-31 
March 2010 for this study. This dataset normally covers 
all NHS and private provider admission spells in any 



Box 1: Predictor variables considered 



Sociodemographic variables 

■ Age at admission 

■ Gender 

■ Socioeconomic deprivation status of place of residence 

Temporal and Geographical Factors 

■ Financial year of admission 

■ Season of admission 

■ Distance of place of residence from hospital 

Health and Health Service Factors 

■ Primary diagnosis at admission 

■ Charlson comorbidity index at admission 

■ Treatment specialty 

■ Patients' GP practice location 




hospital in England under the care of a consultant, 
midwife or nurse. 19 A provider spell is the time that a 
patient stays with one hospital care provider from admis- 
sion to discharge, transfer or death. 20 

We extracted for analyses records of admission spells 
for patients resident in Blackpool aged 40 years or older 
admitted to BTH with a primary diagnosis of COPD 
(ICD-10 code J40^J44) 21 and with no primary proce- 
dures (invasive or non-invasive) carried out or recorded. 
NHS numbers were anonymised and used as unique 
identifiers for admissions attributable to a patient. 

Patients were also grouped into four age bands (40-59, 
60-69, 70-79 and 80+). To identify the subclassifications 
of COPD conditions, ICD-10 codes were linked to 
their respective descriptions. 21 We estimated the 
Charlson comorbidity index 22 using the secondary diag- 
noses recorded at admission. We subsequently grouped 
the patients into three indices bands namely 1, 2 and 3+. 

Postcodes of places of residence were linked to the 
2010 English Index of Multiple Deprivation (IMD), 28 
which are composite scores derived from seven depriv- 
ation domains, using the lower super output codes. The 
scores were grouped into the national IMD quintiles. 
The higher the score, the more deprived an area is. We 
classified the seasons of admission using the Office for 
National Statistics (ONS) classifications used for the 
excess winter deaths calculations. 18 

To identify patient's general practitioner (GP) practice, 
practice codes were linked to their respective locations and 
classified as from Blackpool or from outside Blackpool 
because there were very small numbers of admissions 
recorded for some of the practices. We estimated the short- 
est distance (in kilometres) between patients' places of 
residence and the hospital using the geographical 
co-ordinates (Easting and Northing) associated with post- 
codes. These distances were categorised into quintiles. 



Outcomes 

In this study we sought to identify the factors associated 
with the LOS among Blackpool COPD patients. The 
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factors considered in the analyses are listed in box 1. The 
main outcome measures were the median and mean 
length of hospital stay, and HRs of being discharged 
earlier from hospital associated with the variables. 

Statistical analyses 

We calculated LOS from the date of admission and the 
date of discharge. Discharges to usual places of resi- 
dence, hospices and care homes were considered the 
endpoints for each admission spell while admissions 
resulting in deaths and those resulting in transfers to 
other NHS hospital care providers were censored. 

We used the Kaplan-Meier method to estimate median 
and restricted mean LOS because the LOS data were 
positively skewed 24 and because of the presence of cen- 
soring. 25 To identify factors independendy associated 
with LOS we used Cox proportional hazard model for 
both univariate and multivariate analyses. We tested the 
proportional hazard assumption using the Schoenfeld 
residuals test. 26 The assumption did not hold for gender 



in both univariate and multivariate models (table 1) 
hence we used the stratified Cox model (stratified by 
gender) for the final multivariate model that included 
all the remaining variables. We also adjusted for cluster- 
ing at patient level because there were multiple admis- 
sions recorded for some of the patients over the period. 
The results from the multivariate model containing all 
the main variables were reported as adjusted HRs. 

We added interaction terms for age group and IMD 
quintile (age groupxIMD quintile) and age group and 
season of admission (age groupxseason) to the final 
model in turn to assess any interaction between these 
factors. We also evaluated interactions between season 
and level of deprivation (seasonxIMD quintile) , season 
and distance from hospital (seasonxdistance), and 
patient's age and distance from hospital (agexdistance). 

We used the likelihood ratio test to evaluate the contri- 
bution of the interaction terms to the overall fit of the 
models. We did not use the likelihood ratio test for the 
main effects to identify those that contribute significantly 



Table 1 Variations in length of stay by patient sociodemographic variables 





Number of 


Median 














p For 




admissions 


LOS 


Mean LOS 


Unadjusted 




p For PH 


Adjusted HRf 




PH 


Variable 


(%) 


(95% CI) 


(95% CI) 


HR (95% CI) 


P 


test 


(95% CI) 


P 


test* 


Age group 


















40-59 


410 (17.0%) 


5 (5 to 6) 


6.96 


1 






1 




_ 








(6.13 to 7.80) 














60-69 


741 (30.7%) 


5 (5 to 6) 


7.11 


0.95 


0.6010 


0.5842 


0.90 


0.1830 


0.8206 








(6.57 to 7.65) 


(0.79 to 1.15) 






(0.77 to 1 .05) 






70-79 


788 (32.7%) 


6 (6 to 7) 


10.92 


0.71 


<0.0001 


0.7668 


0.69 


<0.0001 


0.2508 








(9.77 to 12.07) 


(0.58 to 0.86) 






(0.59 to 0.82) 






80+ 


471 (19.5%) 


7 (7 to 8) 


14.51 


0.58 


<0.0001 


0.9745 


0.59 


<0.0001 


0.8858 








(11.84 to 


(0.48 to 0.71) 






(0.50 to 0.70) 












17.18) 














Total 


2410 (100%) 


6 (6 to 6) 


9.79 






0.6999 














(9.09 to 10.50) 














Gender 




















Male 


1152 (47.8%) 


5 (5 to 6) 


9.46 


1 


















(8.25 to 10.66) 
















1258 (52.2%) 


7 (6 to 7) 


10.09 


0.85 


0.0070 


<0.0001* 








Female 






(9.33 to 10.85) 


(0.76 to 0.96) 












Total 


2410 (100%) 


6 (6 to 6) 


9.79 






<0.0001* 














(9.09 to 10.50) 














IMD quintile^ 




















51 (2.1%) 


4 (3 to 5) 


5.29 


1 






1 






Second 






(4.11 to 6.46) 














Third 


311 (12.9%) 


6 (5 to 6) 


10.08 


0.62 


0.0010 


0.2742 


0.65 


0.0080 


0.5542 








(8.22 to 1 1 .95) 


(0.46 to 0.83) 






(0.47 to 0.89) 






Fourth 


642 (26.6%) 


6 (6 to 7) 


11.76 


0.57 


<0.0001 


0.5780 


0.57 


<0.0001 


0.3118 








(9.85 to 13.67) 


(0.44 to 0.73) 






(0.42 to 0.77) 






Fifth 


1406 (58.3%) 


6 (6 to 6) 


8.95 


0.65 


0.0010 


0.9077 


0.59 


0.0010 


0.4996 








(8.30 to 9.60) 


(0.51 to 0.84) 






(0.43 to 0.82) 






Total 


2410 (100.0%) 


6 (6 to 6) 


9.79 (9.09 to 






0.2134 














10.49) 















*p for PH test in the exploratory multivariate model was <0.0001 . 
LOS, length of stay; PH test, proportional hazard test. 
tHR adjusted for all the main variables. 

tThere are no first Index of Multiple Deprivation national quintile areas in Blackpool. 
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to the fit of the model because it was our aim to describe 
any association between these variables and LOS. 

Only the results for the variable involved in the inter- 
action terms were displayed because the results for the 
other main variables were the same as in the final 
model without the interaction terms. We set the statis- 
tical significance level for all the analyses at 5%. We used 
Stata V.10 for all the analyses. 

RESULTS 

Admission outcomes 

There were 2410 admissions meeting the inclusion cri- 
teria over the period. These admissions were attributed 
to 1172 COPD patients, giving an average of 2.06 admis- 
sions per patient. In all 2226 (92.4%) of the admissions 
were discharged to usual places of residence, 12 (0.5%) 
were discharged to local authority residential accommo- 
dation, while 7 (0.3%) were to local authority-run care 
homes and 3 (0.1%) were to local authority-run hos- 
pices. A total of 147 (6.1%) of the admissions resulted 
in deaths while 15 (0.7%) were transferred to other 
NHS hospital care providers. The median LOS was 
6 days (95% CI 6 to 6) while the mean was 9.8 days 
(95% CI 9.1 to 10.5). 

Length of stay and sociodemographic variable 

Table 1 shows the variation in LOS by patient sociode- 
mographic variables. Mean age at admission was 
70.1 years (95% CI 69.7 to 70.5). The minimum age was 
40 and the maximum 100 years. The highest proportion 
of the admissions (32.7%) was for patients aged between 
70 and 79 years while the lowest (17.0%) was for patients 
aged between 40 and 59 years. There were more admis- 
sions for females (52.2%) compared with males. More 
than half (58.3%) of the admissions were for patients 
from the most deprived fifth quintile areas. 

The LOS was significantly longer for those aged 
80+ years compared with those aged 40-59 years. 
Compared with those aged 40-59 years, those aged 80+ 
years were 41% less likely to be discharged earlier from the 
hospital (adjusted HR 0.59; p=<0.0001) while those aged 
70-79 years were 31% less likely (adjusted HR 0.69; 
p=<0.0001). There was no significant difference between 
the likelihood of earlier discharge for those aged 
40-59 years and 60-69 years (adjusted HR 0.90; p=0.1830). 

Females had a slightly longer LOS compared with males 
but the difference was not statistically significant for both 
the mean and median values. The unadjusted HR showed 
that females were 15% less likely to be discharged earlier 
compared to males (unadjusted HR 0.85; p=0.0070). 
However because the proportional hazard assumption did 
not hold for gender (p=<0.0001), the multivariate model 
was stratified by gender. No results for gender were pre- 
sented for the stratified model because stratified Cox 
model does not report HRs for stratified variables. 

There was a significant association between deprivation 
status of place of residence and LOS with those in the 



more deprived areas tending to stay longer on admission. 
Compared with second IMD quintile areas, those from 
fifth quintile areas were 35% less likely to be discharged 
earlier (adjusted HR 0.65; p=0.0010), those from the 
fourth IMD quintile areas were 43% less likely to be dis- 
charged earlier (adjusted HR 0.57; p<0.0001), while those 
from the third quintile areas were 38% less likely 
(adjusted HR 0.62 p=0.0010). 

LOS variations by temporal and geographical factors 

Table 2 shows the variations in LOS by temporal and 
geographical variables. The highest number of admis- 
sions was recorded in 2007/08 financial year while the 
lowest was recorded in 2006/2007. There were more 
admissions during December-March (38.1% of all the 
admissions) compared to the other two seasons. 

The average distance of the places of residence from 
the hospital was 2.99 km (95% CI 2.94 to 3.05) with the 
furthest distance being 6.32 km and nearest 0.32 km. 
Half of the patients were from within 2.75 km (ie, the 
median distance) radius of the hospital. 

There was a statistically non-significant reduction in 
LOS over the period of the study from a median of 
7 days in 2005/2006 to a median of six in 2009/2010. 
The corresponding mean LOS figures were 11.5 and 
9.3 days, respectively. The adjusted HR showed that 
patients were 22% more likely to be discharged earlier 
in 2009/2010 compared with 2005/2006 (adjusted HR 
1.22; p=0.0100). Statistically significant increase in the 
likelihood of earlier discharge was noticed as early as 
2006/2007 (adjusted HR 1.20; p: 0.0130); however, this 
increase was not sustained over the rest of the period. 

Though the number of admissions showed seasonal 
variations, the LOS did not. Overall we did not notice 
statistically significant variations in the likelihood of 
earlier discharge across the three seasons. There were 
also no significant variations in the likelihood of earlier 
discharges based on distance from the hospital. 

LOS by health condition and health service factors 

Table 3 shows the variations in LOS by health condition 
and health service variables. Majority (90.8%) of the 
patients were registered with Blackpool CPs. Just over 
50% of the admissions were attributed to COPD with 
acute exacerbation, unspecified followed by COPD with 
acute lower respiratory infection (40.2%) . The vast major- 
ity of the admissions were managed by general medicine 
specialty (90%) while 7.1% were managed by accident 
and emergency specialty. Admissions in patients with 
Charlson index of 1 were in the majority (75.9%) while 
those with an index of 3+ were in the minority (7.3%) . 

There were no significant variations in LOS based on 
CP practice location. However, the causes of admissions 
were associated with LOS. Compared with admissions 
attributed to COPD with acute exacerbation, unspecified, 
admissions attributed to COPD with acute lower respira- 
tory infection were 19% less likely to be discharged 
earlier (adjusted HR 0.81; p=<0.0001) while there was no 
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Table 2 Variations in length of stay by temporal and geographical factors 
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Table 3 Variations in length of stay by factors associated with patient's health condition and health service factors 
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Median LOS 


Mean LOS 


Unadjusted HR 




p For 


Adjusted HRf 




p For 


Variable 


admissions (%) 


(95% CI) 


(95% CI) 


(95% CI) 


P 


PH test* 


(95% CI) 


P 


PH test* 


Commissioning PCT 




















Blackpool 


2189 (90.8%) 


6 (6 to 6) 


9.65 


1 


— 




1 












(8.96 to 10.33) 
















221 (9.2%) 


5 (5 to 7) 


10.66 


1 


0.9900 


0.0046 


0.97 


(08300 


0 6738 








(7.88 to 13.43) 


(0.81 to 1 .24) 






(0.76 to 1 .24) 






Total 


2410 (100%) 


6 (6 to 6) 


9.79 




- 


0.0046 






- 








(9.09 to 10.50) 














Primary diagnosis 




















COPD with acute 


1210 (50.2%) 


5 (5 to 6) 


8.15 


• 


— 




1 






exacerbation, unspecified 






(7.57 to 8.72) 














COPD with acute lower 


968 (40.2%) 


7 (6 to 7) 


11.29 


0.79 


<0.0001 


000)57 


0.81 


■ 0 0001 


08098 


respiratory infection 






(10.15 to 12.43) 


(0.73 to 0.86) 






(0.75 to 0.88) 






Other COPD 


232 (9.6%) 


5 (4 to 6) 


11.56 


0.91 


0.2790 


0.0047 


0.93 


0.3570 


0.0300 








(7.13 to 15.98) 


(0.76 to 1 .08) 






(0.78 to 1 .09) 






Total 


2410 (100%) 


6 (6 to 6) 


9.79 


— 




0.0163 






— 








(9.09 to 10.50) 














Charlson index 




















1 


1829 (75.9%) 


6 (5 to 6) 


8.93 


1 






1 




- 








(8.33 to 9.53) 














2 


404 (16.8%) 


7 (6 to 7) 


12.22 


0.82 


0 0050 


000)38 


0.87 


0 0390 


0.2730 








(9.58 to 14.86) 


(0.72 to 0.94) 






(0.76 to 0.99) 






3+ 


177 (7.3%) 


8 (7 to 8) 


13.27 


0.73 


<0.0001 


0.7163 


0.81 


0.0080 


0.2212 








(8.76 to 17.78) 


(0.63 to 0.86) 






(0.69 to 0.95) 






Total 


2410 (100%) 


6 (6 to 6) 


9.79 




- 


0.9352 




- 


- 








(9.09 to 10.50) 














Specialty of admission 




















Accident and emergency 


171 (7.1%) 


3 (1 to 5) 


5.87 (4.26 to 


• 


- 


- 


1 


- 


- 








7.48) 














General medicine 


2170 (90%) 


6 (6 to 6) 


10 (9.24 to 10.75) 


0.64 


<0.0001 


0.1253 


0.67 


0.0020 


0.0771 










(0.50 to 0.82) 






(0.52 to 0.86) 






Other 


69 (2.9%) 


6 (5 to 7) 


9.21 


0.64 


0.0060 


0.1085 


0.63 


0.0060 


0.0254 








(6.39 to 12.04) 


(0.47 to 0.88) 






(0.45 to 0.88) 






Total 


2410 (100%) 


6 (6 to 6) 


9.79 






0.2470 














(9.09 to 10.50) 














COPD, chronic obstructive pulmonary disease; LOS, length of stay; PH test, proportional hazard test. 












*HR adjusted for all the main variables. 



















Factors influencing length of hospital stay among COPD patients 



statistically significant difference between the likelihood 
of earlier discharge for admissions due to unspecified 
causes of exacerbations and admissions attributed to 
other COPD problems (adjusted HR0.93; p=0.3570). 

Specialty of admission was also significandy associated 
with LOS with admissions managed by accident and 
emergency specialty being discharged the earliest. 
Compared with patients managed by the accident and 
emergency specialties, those managed by general 
medical specialties were 33% less likely to be discharged 
earlier (adjusted HR 0.67; p=0.0020) while those 
managed by other specialties were 37% less likely to be 
discharged earlier (adjusted HR 0.63; p=0.0020). 

Comorbid conditions in patients were associated with 
LOS. Those with Charlson index of 1 experienced the 
shortest LOS while those with Charlson index of 3+ 
experienced the longest. Compared with patients with 
index 1, those with index 2 were 13% less likely to be 



discharged earlier (adjusted HR 0.87; p=0.0390) while 
those with index 3+ were 19% less likely to be dis- 
charged earlier (adjusted HR 0.81; p=0.0080). 

Interactions 

We noticed statistically significant association for two of 
the interaction terms (ie, IMD quintilexseason and age 
groupxdistance) and presented the results for them only 
(table 4). We did not observe significant association for 
the rest. 

In the second IMD quintile areas, admissions during 
August-November were 94% more likely to be dis- 
charged earlier compared with admissions during April- 
July (adjusted HR 1.94; p=0.0150) but in the fourth and 
fifth deprivation quintile areas the effects were reversed. 
In the fourth IMD quintile areas, admissions during 
August-November were 45% less likely to be discharged 
earlier compared with admissions during April-July 



















p For likelihood 


Interaction term* 


Adjusted HR (95% CI) 


p Value 


ratio test 


beasonxlMU 








[2nd quintile]x[Apr-Jul]t 


1 




0.6943 


[3rd quintile]x [Apr-Jul]f 


0.74 (0.51 to 1 .06) 


0.1040 




[4th quintilejx |Apr-JulJT 


0.66 (0.47 to 0.92) 


0.0140 




[5th quintilejx |Apr-JulJT 


0.74 (0.53 to 1 .04) 


0.0800 




[Aug— Nov]x [2nd quintile] 


H C\A / H H A O O H \ 

1.94 (1.14 to 3.31) 


0.0150 




[Dec— Mar]x [2nd quintile] 


■i 0>1 /<"» CO trt O 0"7\ 

1 .^4 (U.OO tO tL.iil) 


0.4880 




[Aug-Nov]x[3rd quintile] 


0.59 (0.31 to 1.11) 


0.1020 




[Dec-Mar]x[3rd quintile] 


0.86 (0.44 to 1 .67) 


0.6500 




[Aug-Nov]x [4th quintile] 


0.55 (0.31 to 0.97) 


0.0380 




[Dec-Mar]x[4th quintile] 


0.87 (0.46 to 1 .65) 


0.6680 




[Aug-Nov]x[5th quintile] 


0.49 (0.28 to 0.85) 


0.0120 




[Dec- Mar]x[5th quintile] 


0.79 (0.43 to 1 .47) 


0.4560 




Agexdistance 








[40-59]x[Distance-q1]+ 


1.00 




0.0055 


[60-69]x[Distance-q1]f 


0.79 (0.60 to 1 .05) 


0.1050 




[70-79]x[Distance-q1]+ 


0.52 (0.38 to 0.72) 


<0.0001 




[80+]x[Distance-q1]t 


0.44 (0.31 to 0.63) 


<0.0001 




[40-59]x[Distance-q2] 


1 .03 (0.78 to 1 .35) 


0.8570 




[40-59]x[Distance-q3] 


0.95 (0.69 to 1 .29) 


0.7310 




[40-59]x[Distance-q4] 


0.55 (0.39 to 0.76) 


<0.0001 




[40-59]x[Distance-q5] 


0.84 (0.52 to 1 .37) 


0.4930 




[60-69]x[Distance-q2] 


1.19 (0.82 to 1.74) 


0.3620 




[60-69]x[Distance-q3] 


0.96 (0.65 to 1 .43) 


0.8470 




[60-69]x[Distance-q4] 


1.54 (1.03 to 2.31) 


0.0360 




[60-69]x[Distance-q5] 


1.13 (0.64 to 2.00) 


0.6790 




[70-79]x[Distance-q2] 


1.11 (0.74 to 1.69) 


0.6130 




[70-79]x[Distance-q3] 


1 .45 (0.89 to 2.35) 


0.1350 




[70-79]x[Distance-q4] 


1.68 (1.10 to 2.58) 


0.0170 




[70-79]x[Distance-q5] 


1.57 (0.88 to 2.82) 


0.1300 




[80+]x[Distance-q2] 


1.17 (0.75 to 1.84) 


0.4820 




[80+]x[Distance-q3] 


1.13 (0.70 to 1.83) 


0.6160 




[80+]x[Distance-q4] 


2.25 (1.41 to 3.60) 


0.0010 




[80+]x[Distance-q5] 


1.33 (0.76 to 2.35) 


0.3200 




*See http://128.97.141.26/stat/stata/webbooks/logistic/chapter2/default.htm for a detailed explanation of interpreting interaction terms. 


tBaseline risk for within subgroup comparisons. 
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Factors influencing length of hospital stay among COPD patients 



(adjusted HR 0.55; p=0.0380) while in the fifth IMD 
quintile areas, they were 51% less likely to be discharged 
earlier compared with admissions during April-July 
(adjusted HR 0.49; p=0.0120). We did not find signifi- 
cant variations in the likelihood of earlier discharge 
between admissions in April-July and December-March 
across the second, fourth and fifth quintiles. There were 
no statistical significant seasonal variations in the LOS 
for admissions from the third deprivation quintile areas. 

Across all age groups, admissions from the fourth dis- 
tance quintile areas (ie, within 3.1-4.2 km radius of the 
hospital) were significantly more likely to be discharged 
earlier compared with those living within the first quintile 
distance of the hospital. This effect increases with age with 
those aged 80+ and living within fourth distance quintile 
of the hospital being more than two times more likely to 
be discharged earlier compared with their counterparts 
living within the first distance quintile of the hospital 
(adjusted HR 2.25; p=0.0010). The interaction term signifi- 
cantly contributed to the fit of the model ( p=0.0055) . 

DISCUSSION 

Main findings in relation to the literature 

LOS for COPD patients have reduced over the period of 
our study. This was in keeping with national and inter- 
national trends which have also been associated with 
increasing re-admission rates 8 12 27 raising the question 
whether there could be an optimum LOS for COPD 
patients. Though some have suggested 7 days as an 
optimum associated with the lowest re-admission rates 11 
others found no association between re-admission rates 
and LOS. 28 The apparent inconsistencies in the evidence 
about re-admission risk from existing literature may be 
due to developments in postdischarge management. 11 

Older patients were more likely to stay longer on admis- 
sion. This may be because they might have been in poorer 

27 

functional states and these were not accounted for in 
our study. The literature also indicates that LOS tended to 
increase with age. 29 30 Old age is also associated with 
general physical frailty that could delay recovery. 

We found that males were more likely to be dis- 
charged earlier than females though the proportional 
hazard assumption did not hold. There is some evidence 
from the existing literature 11 31 supporting our observa- 
tion though the reasons for this are not apparent. 

We also found that admissions from more deprived 
areas were more likely to stay longer in hospital possibly 
due to them having more severe underlying health con- 
ditions. COPD admissions are most frequent in deprived 
areas. 17 32 The combination of frequent admissions and 
prolonged LOS will accentuate the cost burden of 
COPD in deprived areas. 

Though we did not notice an overall seasonal variation 
in LOS, we observed seasonal variations across the 
deprivation quintiles. Those in the most deprived areas 
were less likely to be discharged early during August- 
November. This may be a reflection of the need to 



involve other services (eg, social services) in the care of 
the patients from more deprived areas in the run up to 
the winter. There is some evidence showing that the 
need for social work intervention may be linked to pro- 
longed LOS in patients with acute exacerbations. 33 

Exacerbations due to infective causes were associated 
with longer LOS in our study. This is supported by find- 
ings in the UK. 34 The risk of respiratory infection is 
highest during the winter. 35 36 Infective causes of exacer- 
bations are therefore expected to be higher in the 
winter and may be associated with longer LOS. It is thus 
surprising that our data did not show any significant sea- 
sonal variation in both univariate and multivariate ana- 
lyses though such variations had been noticed 
elsewhere. 37 Implementation of supported discharge 
schemes in Blackpool and hospital bed management 
pressures during the winter might have contributed to 
the attenuation in seasonal variations. 

Patients treated by general medical specialties were 
most likely to stay longer. This may be due to the 
complex comorbid conditions managed by these special- 
ties. There is however evidence suggesting that LOS are 
influenced by managing physician attributes and the 
quality of care organisation. 30 38 Patients admitted 
under care of the elderly physicians have been found to 
be less likely to enter early discharge schemes. 29 

Comorbidities were independent predictors of LOS in 
our study. This has been found by others as well. 39 40 It 
is reasonable to expect LOS to be longer in patients 
with multiple underlying health conditions some of 
which may be related to COPD. 

Generally the distance of place of residence from hos- 
pital was not associated with LOS but analysis based on 
interaction between age group and distance showed that 
across all the age groups, those living within 3.1-4.2 km 
radius of the hospital were most likely to be discharged 
earlier. This may be a feature associated with service con- 
figurations and possibly the effects of living conditions 
not adequately captured by the IMD scores. Though we 
have not come across any study on the impact of distance 
of patients' places of residence from hospital on COPD 
admission outcomes, Purdy et al 1 found that distance to 
the nearest emergency department was significantly asso- 
ciated with risk of admission for respiratory condition. 

Implications for practice 

Those from more deprived areas may suffer more seriously 
from COPD and comorbidities and may make proportion- 
ately higher demand on services. To reduce these pres- 
sures, improving access to measures aimed at reducing 
exacerbations, for example, pulmonary rehabilitation ser- 
vices, may need to be better focused on these areas. 

Effects of early discharge schemes on seasonal varia- 
tions in LOS, and admission and re-admission frequen- 
cies will need to be evaluated to ensure optimal care is 
provided across all the seasons and that the 'revolving 
door' phenomenon 11 is not operating. 
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LOS is routinely used by health authorities to deter- 
mine efficiency in care provision. Given the fact that no 
optimum LOS has been identified for COPD patients, if 
LOS is used without case-mix considerations it could 
lead to inappropriate inferences about efficiency. 

Limitations of the study 

Incompleteness and inaccuracy are the main problems 
associated with routine data sources. Even though there 
have been concerns about the quality of routine hospital 
data, data quality in the UK has improved greatly over 
the years. 41 The dataset we used is one of the standard 
datasets used for performance monitoring, reconcili- 
ation and payments in the National Health Service and 
supports the implementation of the UK Department of 
Health payment by result policy. 42 

This study was on a patient population in a defined geo- 
graphical area. This has implications for its generalisability 
as factors such as patient profile, seasonal influences and 
service configuration could influence outcomes. 

Other likely predictors of LOS, for example, severity 
of the primary diagnoses and comorbidities, perform- 
ance status and availability and quality of other health 
and social services were not included in the regression 
models because they were not available in the dataset. 
We did not include admission method and discharge 
destination because very few patients were in some of 
the categories to permit appropriate analyses. 
Additionally there was no information on the discharge 
destination (ie, where they would have been discharged 
to should their clinical management achieve the 
outcome of being fit enough to leave hospital) for those 
whose data were censored. We used Charlson index to 
model the effect of comorbidities. Some have observed 
that the exclusive use of the index may underestimate 
comorbidities in COPD patients. 43 

We used the 2010 IMD, a population level measure of 
deprivation, to assess the association between level of 
deprivation and LOS. Though IMD is widely used in 
studies to assess the effect of socioeconomic deprivation 
on health outcomes, it could predispose to ecological 
fallacy because associations at the population level may not 
necessarily represent associations at an individual level. 44 

Finally readers need to be cautious in interpreting 
these associations between LOS and the factors of inter- 
est because the associations do not necessarily imply 
causation. 

CONCLUSION 

LOS among COPD patients have reduced over the 
period of the study. Age, deprivation, Charlson index, 
specialty of admission and cause of exacerbations were 
independently associated with LOS. We did not find sig- 
nificant associations between LOS and season of admis- 
sion and distance of place of residence from hospital 
but there were significant variations in LOS for these 
variables based on selected patient characteristics. 
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